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We explore CP violation effects on the neutrino propagation in dense environments, such as in 
core-collapse supernovae, where the neutrino self-interaction induces non-linear evolution equations. 
We demonstrate that the electron (anti)neutrino fluxes are not sensitive to the CP violating phase 
if the muon and tau neutrinos interact similarly with matter. On the other hand, we numerically 
show that new features arise, because of the non-linearity and the flux dependence of the evolution 
equations, when the muon and tau neutrinos have different fluxes at the neutrinosphere (due to 
loop corrections or of physics beyond the Standard Model) . In particular, the electron (anti)neutrino 

'-^', probabilities and fluxes depend upon the CP violating phase. We also discuss the CP effects induced 

^i: . by radiative corrections to the neutrino refractive index. 

^ ' PACS numbers: 97.60.Bw,14.60.Pq,11.30.Er 
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w ; I. INTRODUCTION 

^^ . One of the major open issues in neutrino physics is the possible existence of CP violation. As with the recent 

,_i' crucial experimental discoveries on neutrino oscillations, the answer to this question has fundamental implications 
,^ , in high-energy physics, astrophysics and cosmology e.g. to understand the matter versus anti-matter asymmetry in 
Oh' the Universe. The observation that the weak interaction violates the CP symmetry in the quark sector was first 
Q [ established in 1964 11|. Future strategics to search for CP violation in the lepton sector depend upon the actual value 
%^ • of one yet unknown neutrino oscillation parameter, i.e. Oi^ [3, and require long term accelerator projects producing 
^ " very intense neutrino beams [3| . It is therefore essential to explore alternative avenues to get clues on this fundamental 
Ci . question, such as indirect effects in dense environments like core-collapse supernovae. 

Core-collapse supernovae emit about 10^^ erg as neutrinos of all flavours during their rapid gravitational collapse. 
^sl ■ Such neutrinos might play a role on the two major supernova unsolved problems, namely understanding how the 
^ ' explosion finally occurs and where the nucleosynthesis of the heavy elements, produced during the r-process, takes 
OO , place. While neutrinos from a massive star were first observed during the SN1987A explosion, future observations 
^~~^ • of (extra)galactic or relic supernova neutrinos will help unravelling supernova physics and/or unknown neutrino 
^" I properties. For example, in [J, |5|, [y, |7| the imprint of the shock wave on the neutrino time signal is investigated; 
; . while avenues for extracting information on the third neutrino mixing angle are discussed in [8J, |9|. This searches 
P^ ' require advances in the modelling of supernova dynamics, of neutrino propagation in dense environments and of our 
^^ , knowledge on neutrinos. 

Impressive developments are currently ongoing in our understanding of neutrino propagation in dense matter. While 

solar experiments |10l . I 111 . Il2l . Il3| have beautifully confirmed the oscillation enhancement induced by the coupling with 

?*■ ' matter - the Mikheev-Wolfenstein-Smirnov or MSW effect [3, [I^ "j recent theoretical investigations have shown that 

k>( [ the inclusion of the neutrino self-interactions in dense environments introducing a non-diagonal neutrino refractive 

'^ ■ index [l^| gives rise to a wealth of new phenomena, as first pointed out in [l3| . Various regimes have been identified 

d ' : the synchronized one [1^, [I^ , the bipolar oscillations [l^, [20| and the spectral split phenomenon [2lj, [231 ■ Since 

numerical calculations become more involved, analytical treatments for the three flavour case are being proposed (see 

e.g. |23l). The importance of the loop corrections to the neutrino refractive index, the Vf^r potential [2J|, is underlined 

in [25[. Moreover constraints on neutrino mixing from shock re-heating and the r-process nucleosynthesis, including 

the neutrino-neutrino interaction, arc investigated in [26|, |27[. The impact of the neutrino-neutrino interaction on the 

electron fraction relevant for the r-proccss is investigated in [2^, [231 ■ 

In a previous work [301 we have investigated the CP effects on the neutrino fluxes, on the electron fraction in a 
supernova (relevant for the r-process) as well as the possible impact on the supernova neutrino signal in an observatory 
on Earth. In particular we have shown analytically that no effects can be found on the electron (anti)neutrino fluxes, 
when muon and tau neutrino have the same fluxes at the neutrinosphere; while significant effects are obtained 
numerically on the fluxes when they differ. The calculations in (30| are obtained considering interaction with matter 
at tree level only. In this work we explore for the first time CP violation effects on the neutrino propagation in 
dense environments including the standard MSW effect, the neutrino self-interactions and the V^r refractive index. 
We first show analytically that if the muon and tau neutrinos interact similarly, the electron (anti) neutrinos are not 
sensitive to the CP violating phase, even in presence of neutrino self-interactions. This result is general and valid for 
any matter density profile and/or initial neutrino luminosity. We present numerical results, obtained within the three 
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flavour formalism, on the neutrino oscillation probabilities and fluxes within the star. In particular we show that 
both the probabilities and the fluxes become sensitive to the CP violating phase if muon and tau neutrino interact 
differently with matter (e.g. because of loop corrections or physics beyond the Standard Model.) The paper is 
structured as follows. Section II presents the theoretical framework for describing the neutrino propagation including 
the coupling with matter as well as the neutrino-neutrino interaction. Section III gives the analytical and numerical 
results. Conclusions are drawn in Section IV. 

II. THEORETICAL FRAMEWORK 

In a dense environment the non-linear coupled neutrino evolution equations with neutrino self-interactions are given 
by ( we follow the formalism of Ref . [311 ) : 

^uV'.^o = [^0 + H„, + H^^]i:^^ (1) 

at - — 

where ipva denote a neutrino created at the neutrinosphere initially in a flavour state a = e, n, t, Hq = UHyacU^ is the 
Hamiltonian describing the vacuum oscillations Hyac = diag{Ei, E2, E^), -Ei=i,2.3 being the energies of the neutrino 
mass eigenstates, and U the unitary Maki-Nakagawa-Sakata-Pontecorvo matrix 

/I \ / ci3 si3e~'^\ / ci2 S12 0\ 

U = TssTigTia = C23 S23 1 -S12 C12 , (2) 

VO -S23 C23/ V-si3e'* ci3 / V 1/ 

Cij = cosOij {sij = sinOij) with ^12, ^23 £^nd ^13 the three neutrino mixing angles. The presence of a Dirac S phase 
in Eq.Q renders U complex and introduces a difference between matter and anti-matter. The U matrix relates the 
mass and the flavour basis 

tpu^ = ^UazllJi. (3) 

i 

The neutrino interaction with matter is taken into account through an effective Hamiltonian which corresponds, at 
tree level, to the diagonal matrix Hm = diag{Vc, 0,0), where the Vc{x) = V2GFNe{x) potential, due to the charged- 
current interaction, depends on the electron density N^{x) (note that the neutral current interaction introduces an 
overall phase only). 

The neutrino self-interaction term is 

H,, = V2Gf J2J2 / -0-= (q') (1 - q • (i')dnc.dq' (4) 

where Gf is the Fermi couphng constant, p = pi^^ {~pt ) is the density matrix for neutrinos (antineutrinos) 

P.^^ C.K IV'.J' '/'.>:. , (5) 

q (q') denotes the momentum of the neutrino of interest (background neutrino) and dn^ is the differential number 
density. In the single-angle approximation, that assumes that the neutrinos are all emitted with the same angle, i.e. 
p(q) ~ p{q)j Eq.lO reduces to 

H.. - ^§^Dir/R^) Y.j[p.Aq')L,M') plSq')L,M')W (6) 

with the geometrical factor 



D{r/R,) = \[l-^l-{^YY (7) 



where the radius of the neutrino sphere is Ri, = 10km, and 



L.Ar, E,) _ y3^^^_^^^^^(^) ^ ^ g^p {E^JT^^ v) ^^^ 

where ^2(77) is the Fermi integral, L° and T^^ are the luminosity and temperature at the neutrinosphere. 
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FIG. 1: Neutrino (left) and antineutrino (right) oscillation probabilities in three flavours, as a function of the distance from 
the neutron-star surface (10 km), including the neutrino-neutrino interaction and V^r refractive index. The different curves 
correspond to electron (anti)neutrinos (dot-dashed), muon (solid) and tan (dashed) (anti)neutrinos. The results are obtained 
solving Eas. (|H7|l numerically for a neutrino energy of 5 MeV as an example. The case of inverted hierarchy and small neutrino 
mixing angle On is shown where the neutrino self-interaction effects are particularly impressive : the regimes of synchronized 
and bipolar oscillations can be recognised in the first 100 km. In the case of the electron neutrinos (left figure), the spectral 
split is also apparent. 



III. CP EFFECTS IN PRESENCE OF THE u-u INTERACTION AND V^r REFRACTIVE INDEX 

A. Analytical results 

One way to study under which conditions the electron (anti)ncutrino survival probabihties depend upon the CP 
violating phase S, is to demonstrate that the S dependence of the total Hamiltonian Ht — Hq + H^ + H^i, factorises 
as follows (soil: 



in the T23 basis which is: 



HriS) = SHriS = 0)5*^ 



Ipf, = COs923'lptj.-Sm923tpT, 
'tpr = sin 6*23 -0^1 + COS 6*23 ■0T- 



(9) 
(10) 



Here the whole dependence is in the unitary diagonal matrix S — diag{l, l,e ). In fact, it is straightforward to 
show that for any such Hamiltonian the corresponding evolution operator also factorizes as 



U{S) = 5^t/((5 = 0)5* ^ Hid) = S^H{S = 0)5 

As a consequence one can demonstrate that the electron (anti)neutrino survival probabilities satisfy P(i^e 
0) = P(i^e ^ i^ei <5 = 0) a-iid that the appearance probabilities satisfy : 



P{, 



l^e, S^0)+ P{Vr ^Ve,5^Q)= P{Vt. ^V„5^0)+ P{Vr ^V,,5^ 0). 



Since the v^ {ve) fluxes are given by 



Li,^P{Vr -^ Ve) 



(11) 

(12) 
(13) 



where L,^^ are the neutrino fluxes at the ncutrinosphere, from Eqs. (|12m3p one can see that (j),^^ {(f'Pa) ^^^ '^ot sensitive 
to the CP violating phase if muon and tau neutrinos interact with matter in the same way (i.e. L^ ^ = Lu^). This is 

demonstrated in [30| in the case of the standard MSW case only, i.e. H ~ Hq + Hm- Note that a similar conclusion 
is drawn in [34I, using a different procedure. Besides, from Eqs (|12ll3p . one can see that if the condition Li, ^ L^^ 
is relaxed then (j)y^ and (j)i:^ become dependent on delta, as flrst pointed out in [30|. 
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FIG. 2: Neutrino (left) and anti-neutrino (right) spectra, at 200 km from the neutron-star surface. The different curves 
correspond to : the original Fermi-Dirac distributions for v^ (dotted) and i/^ (solid); the v^ (dashed) and fp (dot-dashed) fluxes 
after the evolution in the star with the neutrino self-interaction. The results are obtained for an inverted hierarchy and a small 
third neutrino mixing angle. While neutrinos show a spectral split, anti-neutrinos undergo full flavour conversion. 



We now show that Eq. ([9]) is indeed satisfied for the total Hamiltonian of Eq. ([T]) including the non- linear H^v term 
of Eq.Q. Let us start with the Liouvillc-Von Neumann equation for the density matrix {h~V): 



i - = [UH.acU^ + H,n + H,,{S), p,„ (<5)], 
at - 



(14) 



To prove our result (that the CP-violating phase can be factorized out of the total Hamiltonian which includes H^,^), 
one has to rotate in the T23 basis, since the S matrix contained in T13 (which can be re written as Tia = S^T^^S ) 
does not commute with T23. We then obtain : 



where 



dSpu„(S)S^ „ + „ t - 4. 4. 






(15) 



(16) 



Let us now consider the evolution equation of the linear combination "^^ Ly^Sp^^ (q, 5)S'^ at a given momentum 
q. At the initial time, this quantity reads, in the r23 basis of Eq. fTU)) . as : 



Y,L,^Sp,M,5,t = 0)S^ = 



(L. 











\ 



(17) 



c\^^^^^ s\^L^^_ C23S23e '\L^^- L^^] 

\ C23S23e*'^(i^^-L^^) sls^i^^ + claP^^ J 

One immediately sees that this quantity does not depend on S if and only if L^, = L^^ . Moreover, the total 

Hamiltonian of Eq. p^ is independent of S at initial time since T^;^Ti2HyacTi2Ti:} + Hm does not depend on 5 (at 
any time) and 



SH,,{t = Q,5)S^ = V2GfY.J{1 -ei-q')[Sp,.At = 0,q')S^ L^Jq') - Sp^Jt = 0,q')S^Lj;Jq')]dq' (18) 



is equal to H^^{t = 0,5 = 0) initially when L^ ^ ~ Ly^(and Li^ ^ = Ljj^). In that case, one can see by recurrence from 
the Liouville-Von Neumann equation Eq. p4|) that the evolution of the term Y^^ L„^Sp„^ (q, 6)S^ is exactly the same 
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FIG. 3: Ratios of the electron neutrino oscillation probabilities for a CP violating phase S = 45° (dotted), 90° (dashed), 135° 
(dot-dashed), 180° (solid) over 5 = 0°, as a function of the distance from the neutron star surface. The left figures shows the 
ratios up to 2000 km while the right figure presents the region between 50 to 100 km where collective effects induced by the 
neutrino self-interaction are maximal. The results correspond to the case of inverted hierarchy and small third neutrino mixing 
angle, for a neutrino energy of 5 MeV. 

as the term ^^ L^^p^^(q,S — 0), since they have the same initial conditions (for any q) and the same evolution 
equations. Indeed, the exact same relation applying at the same time for the anti-neutrino case (where the sign of S 
has to be changed), one simultaneously obtains that at any time: 



hence 



H,,{6} = SH,,{6 = 0)Sl 



HTi6)^SHTiS = 0)S''. 



(19) 



(20) 



Note that the derivation holds both for the multi-angle case Eq.(|3]) and the single-angle case Eq.(IS]). 

This implies that the (anti-)electron neutrino survival probability is independent of S and Ea. ((T2)) is valid, therefore 
01/,, (^) = 4'u^{S = 0) and (j>v^{S) = <i>vS^ = 0), even considering the presence of neutrino-neutrino interaction, if muon 
and tau neutrino fluxes at the neutrinosphere are equal. When the fluxes Li,^ and L^^ are different, the derivation 

does not hold anymore, since ^^^ Ly^Spy^{(\^5)S'^ initially depend on 5. 

In the case where radiative corrections to the neutrino scattering are considered, there is an extra term in Hm = 
diag(Vc, 0, y^r)- For the corresponding Hamiltonian the 5 dependence cannot be factorized anymore and in general 
nothing prevents the electron (anti)ncutrino survival probabilities to be sensitive to the CP violating phase^. These 
cases will be studied in the following. 

B. Numerical results 

The main goal of this section is to investigate numerically the effects that can arise when the factorization Eq.([S]) 
is not satisfied. In particular this occurs in three cases : (i) if the MSW Hq + Hm Hamiltonian does not satisfy 
Eq. pT|) because the V^r refractive index is included; (ii) if the initial conditions Eq. (fT7|) of the H,^^ Hamiltonian are 
not independent on 5, that happens when L^^^ ^ L,j^ , which implies that the neutrino self-interaction term does not 
follow Eg. p^ : (iii) when (i) and (ii) occur. Therefore the results that we present here correspond to the following 
possibilities : 



^ Note that the possible inclusion of nonstandard neutrino interactions in the flavour neutrino mixing as e.g. in |33l| implies that Eq. (9) 
does not hold anymore. 
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FIG. 4: Ratios of the v^ fluxes for a CP violating phase & — 180° over J = 0° as a function of neutrino energy, at 1000 km 
within the star. The curves correspond to the following cases : _ff^^ = and V^r = (dotted), _ff^y j^ and V,it = (dashed) 
and H„^ / and V^^r 7^ (solid). These are obtained with L^^ / L,,^ e.g. T^,^ — 1.05 T,^^. The case H,y„ / and V^t / 
(dot-dashed) with L^^^ = i^^ is also shown. The results correspond to an inverted hierarchy and a small &13. 



a) Hyi, 7^ and V^r 7^ with the condition L^ , = L^^ ; 

b) Hui, ^ and V^t = with the condition L,^ ^ L^,^ ; 

c-} Hvv 7^ and V^t- 7^ with the condition L^ 7^ L^,^ . 

It is important to note that a new feature arises in the a), b) and c) cases : the P(i'e -^ v^) becomes dependent 
on 5. (It is a well known fact that the electron survival probability does not depend on CP violation in vacuum 
and in presence of the interaction of matter at tree level.) In the following we investigate the CP effects by varying 
(5 G [0°, 180°] while we show results for 5 = 180° when at such value the effects are maximal. 

The numerical results we present are obtained by solving the three flavour evolution equation of Eq.(IT]) with a 
supernova density profile having a 1/r^ behavior that fits the numerical simulations shown in [29|. Note that with 
such density profile the region of the first 100 km, where the neutrino self-interaction dominates, is well separated 
from the one of the MSW (high and low) resonances, produced by the interaction with ordinary matter. 

The oscillation parameters arc fixed at the present best fit values [SJ], namely ATO22 = 8 x 10~'^eV^, sin^20i2 — 
0.83 and Ato|3 = 3 x 10^'^ eV^, sin^26'23 = 1 for the solar and atmospheric differences of the mass squares and 
mixings, respectively. For the third still unknown neutrino mixing angle ^13, we take either the present upper limit 
sin^20i3 = 0.19 at 90 % C.L. (L) or a very small value of sin^20i3 = 3 x 10"'' (5) that might be attained at the future 
(third generation) long-baseline experiments [3| . To include the neutrino- neutrino interaction we use the single-angle 
approximation of Eqs.(IB][7]) with Ri, = 10 km, considering that all neutrinos are emitted radially. The neutrino fluxes 
at the neutrinosphere Li,^ Eq.® are taken as Fermi-Dirac distributions with typical average energies of < E^^ > = 
10 MeV, < Ep^ > = 15 MeV and < E,y^ >= 24 MeV with i/^ = t'/^ , ^^r , i^/j , i^r , unless stated otherwise (the chemical 
potentials are assumed to be zero for simplicity). We take the neutrino luminosity L° = lO'^"' erg • s~^. 

Our numerical results in three flavours present the collective oscillations induced by the neutrino-neutrino interac- 
tion, already discussed in the literature (see e.g. [ij, [l^, E, [20, [2l|, [22, [2^, |3al ) . Figure [T] presents the (anti)neutrino 
oscillation probabilities within the star. One recognizes the synchronized regime in the first 50 km outside the neu- 
trinosphere Ri, (assumed here to be equal to the neutron-star surface). In this regime the strong neutrino-neutrino 
interaction makes neutrinos of all energies oscillate with the same frequency so that flavour conversion is frozen, as 
discussed e.g. in [l8|, [36[. When the neutrino self-interaction term becomes smaller, the ordinary matter term starts 
to dominate producing large bipolar oscillations (between 50 and 80 km) that produce strong flavour conversion for 
both neutrinos and anti-neutrinos, in particular for the case of inverted hierarchy, independently of the 613 value 
J20| . Finally neutrinos show complete (no) flavour conversion for energies larger (smaller) than a characteristic energy 
Ec = 7.4 MeV, due to lepton number conservation [21[. This is known as the spectral split phenomenon (apparent 
around 150 km on Figure [U left). (Note that the effect of the partially non-adiabatic MSW low resonance can be 
seen at around 270 km of Figure [Tl right.) 

The neutrino-neutrino interaction might have an important impact on the neutrino spectra as well. If in the case of 
normal hierarchy the flavor evolution of both electron neutrinos and anti-neutrinos are essentially the same as in the 



case where matter only is included, for the case of inverted hierarchy, important modifications are found compared to 
the MSW case [33|- While electron neutrinos swap their spectra with muon and tau neutrinos (Figure [H); the electron 
anti- neutrinos show a complete spectral swapping (Figure [21) • Such behaviours are found for both large and small 
values of the third neutrino mixing angle, in constrast with the standard MSW effect. (Note that anti-neutrinos of 
energies less than 2 MeV have already undergone the MSW low resonance at 200 km, as can be seen from Figure [2j 
right.) 

Let us now discuss the CP violation effects in presence of the neutrino-neutrino interaction-^ and of the loop 
corrections to the neutrino refractive index, with the condition that the muon and tau fluxes at the neutrinosphcre 
are equal {L^ — L^^). Figure [3] shows the ratios of the electron neutrino oscillation probabilities for different 5 values, 
as a function of the distance within the star. A 5 MeV neutrino is taken, as an example. One can see that the 5 effects 
are at the level of 1 %. Note that the presence of H^^ with VJ^,- amplifies these effects that are at the level of less than 
0.1% and smaller, when V^r only is included"^. One can also see that in the synchronized regime the CP effects are 
"frozen" while they develop with the bipolar oscillations. Similar modifications are also found in the case of electron 
anti- neutrinos, with effects up to 10% for low energies (less than 10 MeV). Note that the latter might be partially 
modified in a multi-angle calculation, since it has been shown that the decoherence effects introduced by multi-angles 
modify the electron anti- neutrino energy spectra in particular at low energies [37|. Multi-angle decoherence is also 
discussed in [3l|, [3a, Sll, 143 • To predict how these effects modify the numerical results presented in this paper would 
require a full multi-angle calculation. 

The modifications induced by 5 on the electron neutrino fluxes are shown in Figure [H for the aj, 6J, and c) cases, in 
comparison with a calculation within the MSW effect at tree level only as investigated in a previous work [30| ■ Figure 
[S] shows how the CP effects evolve as a function of the distance from the neutron-star surface for the a) and c) cases. 
To differentiate the muon and tau neutrino fluxes at the neutrinosphcre here we take as an example T^, = 1.05 T^^ 
(note that in [301 differences of 10% are considered). In general, we have found that the inclusion of the neutrino 
self-interaction in the propagation reduces possible effects from 5 compared to the case without neutrino-neutrino 
interaction, as can be seen in Figure [H In all studied cases both for v^ and i>e we find effects up to a few percent at 
low neutrino energies, and at the level of 0.1% at high energies (60 - 120 MeV). Our numerical results show deviation 
at low energies that can sometimes be larger than in absence of neutrino self-interaction (Figure [5]), those at high 
energies turn out to be much smaller. This effect of the neutrino-neutrino interaction might be due to the presence 
of the synchronized regime that freezes possible flavour conversion at initial times and therefore also reduces the 
modifications coming from a non zero CP violating phase at later times. 

IV. CONCLUSIONS 

We have investigated possible effects of the CP violating phase on the neutrino propagation in dense matter when 
interaction with matter without/with loop corrections and the neutrino self-interaction are included. Our analytical 
results demonstrate that, at tree level, even when the neutrino- neutrino interaction is included there are no CP 
violating effects on the electron (anti)neutrino fiuxes in the star unless u^ and v-r fluxes differ at the neutrinosphcre. 
If such condition is not satisfied, a totally new feature arise, namely that the electron (anti)neutrino oscillation 
probabilities (and fiuxes) become sensitive to the CP violating phase 5. The latter is also true when the loop 
corrections to the refractive index are included. We find numerically that, in most cases studied, the modifications 
introduced by the CP violating phase are larger (smaller) at low (high) energies than in the case where the neutrino- 
neutrino interaction is not included, and at the level of a few percent. We also find numerically that, even assuming 
that the muon and tau neutrinos have the same fiuxes at the neutrinosphcre, the CP effects induced by V^t only are 
amplified by the neutrino self-interactions up to several percent. 

The authors acknowledge the support from "Non standard neutrino properties and their impact in astrophysics 
and cosmology" , Project No. ANR-05-JC JC-0023. 



^ Note that a comment is made in 1381 . 13911 on the 5 effects on the neutrino fluxes in the presence of the neutrino self-interaction in a 

core-coUapse supernova. 
^ Note that it was first pointed out in [4011 that the inclusion of the V^t refractive index renders the electron neutrino survival probability 

slightly 6 dependent. 
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FIG. 5: Ratios of the v^ fluxes for a CP violating phase 5 = 180° over 5 = 0° as a function of neutrino energy. They correspond 
to inverted hierarchy and small ^13 and different distances from the neutron star surface, i.e. 200 km (dotted), 500 km (dashed), 
750 (dot-dashed), 1000 (solid). The results include the v-v interaction and the V^,- refractive index. They are obtained using 
equal v^^ and Vr fluxes at the neutrinosphere (left) or taking Ty^ = 1.05 T^^ (right). For the v^ fluxes deviations up to 10% are 
found at energies lower than 10 MeV. 
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